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<N ■ ABSTRACT 

We report the results of a series of three-dimensional (3-D) simulations of 
the deflagration phase of the gravitationally confined detonation mechanism for 
Type la supernovae. In this mechanism, ignition occurs at one or several off- 
center points, resulting in a burning bubble of hot ash that rises rapidly, breaks 
! through the surface of the star, and collides at a point opposite breakout on 

| the stellar surface. We find that detonation conditions are robustly reached in 

our 3-D simulations for a range of initial conditions and resolutions. Detonation 
t — | conditions are achieved as the result of an inwardly-directed jet that is produced 

by the compression of unburnt surface material when the surface flow collides 
with itself. A high-velocity outwardly-directed jet is also produced. The initial 
conditions explored in this paper lead to conditions at detonation that can be 
expected to produce large amounts of 56 Ni and small amounts of intermediate 
mass elements. These particular simulations are therefore relevant only to high 
luminosity Type la supernovae. Recent observations of Type la supernovae imply 
a compositional structure that is qualitatively consistent with that expected from 



these simulations. 
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Introduction 



Type la supernovae have received increased interest because of their importance as 
"standard candles" for cosmology. Observations using Type la supernovae as standard 
candles have revealed that the expan sion rate of the universe is accelerating and have 
led to the discovery of "dark energy" (IRiess et al.l Il998l ; iPerlmutter et al.lll998l ). But the 
way in which Type la supernovae explode is not completely understood. The current 



leading paradigms for the explosion mechanism are ( 1) pure deflagration (IReinecke et al. 



2002b 



Gamezo et al. 



sition (DDT) flKhokhlov 



2003 



Rqpke fc Hillebrandtl l2005h. (2) deflagration to detonation tran 



1991 



pulsational detonation (PD) ( Khokhlo 



Townslev et al. 


2007 


2002; 


Badenes et al. 



1991 



Gamezo Khokhlov fc Oranl 12004: iGamezo et al.l l2005T). (3) 



Bravo et al. 



confined detonation (GCD) ( Plewa. Calder. Sz Lambl 



2004 



20061). and (4) g r avitationally 



Livne et al. 2005; Plewa 



Th ere is increasing ey i dence that a detonation is n eeded (jHoflich et al. 



20061 : IWang et al.ll200d 120071 : iGerardv et al.ll2007l ). as is posited in the 



2007 



last three models. 

A fundamental question has been how the transition to a detonation occurs in a white 
dwarf star [see, e.g., iNiemeyerl (119991 )] . While the DDT, PD, and GCD paradigms incor- 
porate a detonation, all existing DDT simulations invoke the transition to a detonation 
in an ad-hoc fashion and the PD mechanism remains largely unexplored by detailed sim- 
ulations. In contrast, extensive two-dimensional (2-D) cylindrical simulations have shown 
that deton ation conditions are robustly reached in the GCD model for a variety of initial 
conditions rtPlewa. Calder. fc Lambl booi |piewall2007l : iRopke. Wooslev. fc Hillebrandtl 120071 ; 
Townslev et al.ll2007l ). Thus, to date, the GCD mechanism is the only proposed mechanism 
for which a detonation has been demonstrated to arise naturally. 

Howe yer, the achievement of detonat i on co nditions has not been demonstrated in 3-D 
[see, e.g., IRopke. Wooslev. fc Hillebrandtl (120071 )]. This is a concern since the behavior of 
turbulence is different in 3-D than in 2-D, and the cylindrical symmetry of the 2-D simula- 
tions might enhance the focusing of the surface flow that triggers the detonation in the GCD 
model. Hence, a major question that we address in this paper is whether it is possible to 
achieve detonation conditions in a fully 3-D simulati on of the GCD model. In additi o n, the 
large amount of nuclear burning that occurs in the IRopke. Wooslev. &; Hillebrandtl (120071 ) 
simulations appears to play a role in the failure of these simulations to achieve detonation 
conditions in 3-D. In this paper, we report the results of seven 3-D simulations of the GCD 
model for several different sets of initial conditions and two resolutions. We find that the 
conditions for detonation are robustly achieved for these initial conditions. Thus the simu- 
lations reported in this paper address the first point above and provide a counter example 
to the second. 
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The organization of the paper is as follows. We describe the simulation setup in §2 
and the results of the simulations in §3. We discuss the properties of Type la supernovae 
expected in the GCD model and compare the results of our simulations with earlier work in 
§4. Finally, we state our conclusions in §5. 



2. Simulation Setup 



We perform our 3-D simulations of the deflag ration phase of the GCD model usin g 
FLASH 3.0, an adaptive-mesh hydrodynamics code (IFryxell et al.ll2000l ; ICalder et al.ll2002l ). 
The general simulation setup is id entical to that in the 2-D simulations of the GCD model 
described in iTownsley et al.l (120071 ). In particular, the initial model is a 1.38 M & WD with 
a uniform composition of equal parts by mass of 12 C and 16 O. It has a central density of 
2.2 x 10 9 g cm -3 , a uniform temperature of 4 x 10 7 K, and a radius of approximately 2,000 
km. 

The physical thickness of the carbon burning stage of the nuclear flame front in the 
deflagration phase of Type la supernovae is 10~ 4 to 10 3 cm for the densities of interest, and 
is therefore unresolvable in any whole-star simulation. Consequently, a method must be 
used to determine the location and speed of the flame front. Two fairly different methods of 
flame-front tracking have been used in recent studies of WD deflagration. One is the level 
set technique, in which the location of the flame front is calculated based upon the value 
of a smooth field defined on the gr id and propagated w i th an advection equation acting 



in add i tion to the hydrodynamics (IReinecke et al.l Il999l ; IRopke. Niemeyer. fc Hillebrandt 
2006al ; IRopke fc Hillebrandtl 120051 ). Thi s method has been used to study the effect of 



turbulence on the nuclear burning rate (ISchmidt. Niemever. &: Hillebrandtl l2006al) and in 



many s imulations of WD deflagration [see, e.g., 



Reinecke et al 



Jl999h : IRopke fc Hillebrandtl 



( 120051 ); ISchmidt. Niemeyer. fc Hillebrandtl (l2006bl )]. The other method artificially broadens 
the flame front using a reaction pr ogress variable and propagates it using an adyection - 
diffusion-reaction (ADR) equation (jKhokhlovl Il995l ; IVladimirova. Weirs. &: Ryzhikl 120061 ). 



This ADR flame model has been used t o study t 



instability on a propagating flame front (IKhokhlov 



re effect of the Rayleigh - Taylor (R-T) 



previous previous simulations of WD deflagration ( 



Plewa. Calder. k Lambl 12004 ; |Plewall2007l : iTownslev et al.l 120071 ). 



Gamezo et al 



199 5!; IZhang et al.ll2007D and in many 



2003 



Calder et al. 2004 



In this paper, we follow the nuclear flame using a new version (lAsida et al.l 120081 ) of 
ADR flame model. The new prescription uses the Kolmogorov-Petrovski-Piskunov (KPP) 
form of the reaction term in which this term is slightl y truncated, as opposed to the top- 
hat form used previously by ourselves and others [e.g., IKhokhlov! (119951 )]. This new version 
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is nu merically quieter, more stable, and exhibits far smaller curvature effects (lAsida et al. 
20081 ). The thickness of the flame is ~ 4 grid points; more details and the explicit value s 
of the parameters in the ADR flame model tha t we use are given in Townsley et al.l (120071 ). 
We also use a new, acoust ically-quiet version (ITownsley et al.l 120071 ) of the nuclear energy 
release method described in lCalder et al.l (120071 ) that accounts more accurately for the nuclear 
energy released in the flame and in the evolution of nuclear statistical equilibrium (NSE) as 
conditions change within the bubble of hot ash. 

We do not include nuclear burning outside the flame in the simulations reported in thi s 
paper. This approach is the same as that adopted by lRopke. Woosley. fc Hillebrandtl (120 07). 
Were we to have included such burning, both the 2-D simulations reported in lTownsley et al. 
(120071 ) and the 3-D simulations reported in this paper would detonate at th e first instant 
at which detonation condi tions are reached, as we will report in later papers (jMeakin et al. 
2008al ; I Jordan et al.l 120081 ) . Thus including nuclear burning outside the flame would not 
have allowed us to demonstrate that the simulations reported in this paper robustly reach 
conservative conditions for detonation; i.e., that they exceed the temperature and the density 
needed to detonate for a significant period of time. 



We treat the effect of R-T-driven turbulence in the same manner as in ITownsley et al. 

a^l AgrrifAx, where Ax is the 
0.06 is a constant that 



(120071 ); namely, we impose a minimum flame speed s m i n 
grid size of the sim ulation, a = 0.5 is a geomet rical factor, and rrif 



we have calibrated (ITownsley et al.ll2007l . 120081 1 . This is a conservative approach that allows 
the simulation to treat the effects of R-T driven turbulence on resolved scales while ensuring 
that turbulence on the scale of the (artificially broadened) flame thickness does not disrupt 
the flame. We disable this prescription in the truncated cone encompassing the region where 
the surface flow collides with itself to ensure that no unrealistic heating occurs. 

If the nuclear burning rate, as well as the overall dynamics of R-T-driven turbulent nu- 
clear burning, are determined by the behavior at large scales, turbulence at small scales does 
not increase the burning rate and the moderate resolution possible in current simulati ons is 
adequate. That this may be the case is suggested by the results of IZhang et al.l (120071 ). who 
find that the time-averaged rate of buoyancy-driven nuclear burning did not vary when the 
resolution was varied by a factor of four. Because we rely on the resolution of our simulations 
to describe turbulent nuclear burning, and therefore the rate of nuclear burning, we have paid 
close attention to how the rate varies with resolution. We have found no evidence for more 
than a modest variation of the nuclear burning rate due to unresolved behavior - as we discuss 
below. H owever, the appropriat e way to treat turbule n t nuclear burning is an open ques- 
tion |cf . iReinecke et al.l (|2002af) ; IZingale et al.l (120051 ); ISchmidt. Niemeyer. fc Hillebrandt 



( I2006al ); IZhang et al.l ( 20071 )]. and further studies are needed in order to answer definitively 
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this question. 

The core of the star is though t to be convective and therefore to have velocities on 
the l argest scales of ~ 100 km s^ 1 (jWoosley. Wunsch. fc Kuhlenl 12004 Wunsch fc Wooslev 
2004), which are comparable to the laminar flame speed in this region (ITimmes fc Wooslev 
19921 ). The temperature fluctuations in the convective region may lead to one, to a few, or 
to many ignition points. The convective motions are likely to distort the ignition region(s), 
seeding later Rayleigh- Taylor instability modes. However, t he motions are not strong enough 
to destroy the ignition region(s), once born. In addition, iLivne et al.l (120051 ) find that the 
general outcome of off-center ignitions is not strongly affected by the presence of a convective 
velocity field. Finally, we would like to understand the behavior of the simpler case of a single, 
spherical ignition region and zero velocity in the core of the star before considering more 
complicated cases. For all of these reasons, we adopt a single ignition point at a range of 
offset distances from the center of the star as the initial conditions in this paper. We will 
investigate ignition at multiple points and the effects of convective motions in the core of 
the star in future papers. 

We model the ignition region as a spherical bubble of hot ash initially at rest, character- 
ized by an initial radius, rabble , and an initial distance, r fj set , from its center to the center of 
the star along the z-axis. The edge o f the burned reg i on for ms a smooth transition from fuel 
to ash that is ~ 4 zones in width [see iTownsley et al.l (120071 ) for more details] . The density of 
the hot ash is chosen to maintain pressure equilibrium with the surrounding material. The 
surface of the bubble of hot ash correspon ds to the = 0-5 iso surface, where 0i is the flame 
progress variable for 12 C + 12 C burning (ITownsley et al.l 120071 ) . Thus the radius, /"bubble, is 
approximately the radius of the 4>i — 0.5 isosurface. 

Initially, the spherical bubble of hot ash rises slowly, due to its small size and the small 
value of the acceleration of gravity g near the center of the star. The growth of the bubble 



' Vladimirova 


2007; 


Fisher et al. 


2008) 



minimum wavelength for the unstable Rayleigh- Taylor growth of flame surface perturbations, 
s is the laminar flame speed, A = (pf ue i — Pash)/(Pfuei + Pash) is the Atwood number, and g is 
the acceleration of gravity at r g- set . Otherwise, the initial bubble is immediately unstable to 
the growth of Rayleigh- Taylor modes. This is inconsistent with the assumption of a small, 
spherical ignition region. Thus, rabble < A c for all of the simulations reported in this paper , 



except for two we conducted to compare with those of lRopke. Wooslev. fc Hillebrandtl (120071 ) . 



The adaptive mesh refinement criteria we use are chosen to capture the relevant physical 
features of the burning and the flow at reaso nable computat i on ex pense. The criteria for 
refinement are the same as those described in ITownsley et al.l (120071 ). except that we maxi- 
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mally refine a sphere of radius r re fi n e = 1000 km at the center of the star and a truncated cone 
encompassing the region where the flow of hot bubble material over the surface of the star 
collides with itself. The truncated cone has an opening half angle 9 cone = 30°, and extends 
from a radius r cone min = 1500 km to a radius r cone max = 3000 km. 



Results 



We performed a suite of seven 3-D simulations of the deflagration phase of the GCD 
model for initially stationary, spherical flame bubbles, varying the initial bubble radius 
r bubb i c , the initial bubble offset r offset , and the finest resolution. Table 1 lists the initial con- 
ditions and several properties of these simulations. The simulations are denoted by initial 
bubble radius, offset distance, and finest resolution. Thus 18b42o6r denotes a 3-D simulation 
in which the initial bubble radius rbubbie = 18 km, the offset distance r g- sct = 42 km, and 
the finest resolution is 6 km. 

Our 3-D simulations of the GCD model progress similarly to previo us 2-D GCD sim- 



ulations (IPlewa. Calder. fc Lambl 12004 ; iPlewal 120071 ; iTownsley et al. 



20071 ). passing through 



several distinct stages. Initially, the spherical bubble of hot ash that we adopt as our model 
of the ignition region grows at a rate dictated by the laminar flame speed. At ~ 0.2 — 0.3 s of 
simulation time, the radius r bubble of the bubble exceeds A c , the minimum wavelength for the 



unsta ble Rayleigh- Taylor growth of flame surface perturbations (lKhokhlovill995l ; IZhang et al. 



20071 ). When this happens, the top surface of the bubble develops a bulge, and the bub- 
ble quickly evolves into a mushroom-like shape (jCalder et al" 2004 ; Plewa. Calder. fc Lamb 
2004 ; IVladimiroval 120071 ; IPlewal 120071 ; ITownsley et al.l 120071 ). Subsequently, the shape of the 
bubble becomes ever more complex as the critical wavelength A c becomes smaller (as the 
bubble rises and g at the position of the bubble increases) and additional generations of 
smaller features appear as a result of the Rayleigh- Taylor instability. During this time, the 
rate at which the bubble rises increases, and the bubble breaks through the stellar surface 
at ~ 0.8 — 1.2 s. The hot ash in the bubble (which was produced at a range of densities and 
so consists of both iron-peak and intermediate-mass elements) then spreads rapidly over the 
surface of the star, pushing unburned material in the outermost layer of the star ahead of it. 

The mass of the hot ash from the bubble that sweeps over the surface of the star ranges 
from 0.038-0.010 M for offset distances of the initial bubble ranging from 20 km to 100 
km. Its composition ranges from roughly 0.009-0.002 M & of intermediate mass elements and 
roughly 0.007-0.03 M of Fe-peak elements for these offset distances. However, we caution 
that the present simulations do not seek to accurately treat mixing between the ash flow 
across the surface of the star and the underlying outermost layers of the white dwarf, which 
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will significantly affect the final composition of these layers. 



At ~ 1.8 — 2.2 sec, the flow of hot ash collides with itself at the opposite point on the stel- 
lar surface from the place where t he bubble broke out, compresses the unburnt surface layers 



there , and initiates a detonation (IPlewa. Calder. fc Lambl 12004 ; |Plewall2007l ; iTownsley et al. 



20071 ). Figures [T] and [2] illustrate these different stages in the 18b42o6r simulation. The 
image in Figure [3] is identical to that in Figure EJ^d), except we have rotated the star to show 
the very hot (T > 3 x 10 9 K) region at the head of the inward jet that has reached densities 
greater than 2 x 10 7 g cm" 3 , co nditions that exceed conservative criteria for initiation of a 



detonation (INiemeyer &: Woosleylll997l ; iRopke. Woosley. k, Hillebrandtl 120071 ) 



We compare the nuclear energy, E nuc , released during the deflagration phase as a func- 
tion of time for the 18b42o6r and 16b40o8r simulations, and the 25bl00o6r and 25bl00o8r 
simulations, in the left panel of Figure 0} In both cases, the curves for the two different res- 
olutions are in close agreement. This result is expected if buoyancy-driven nuclear burning 
depends mostly on the fluid dynamica l behavior at larger scales. As remarked earlier, such 
a picture is suggested by the results of IZhang et al.l (120071 ). who find that the time- averaged 
rate of buoyancy- driven nuclear burning did not vary when the resolution was varied by a 
factor of four. However, this is an open question, as discussed in §2, and further studies are 
needed in order to answer it. 

We compare the nuclear energy released, E nuc , as a function of time for the 16b20o8r, 
16b40o8r, 16bl00o8r, 25bl00o8r, and 50bl00o8r simulations in the right panel of Figure HI 
The first three form a sequence in which both the time at which the curves flatten and the 
amount of nuclear energy E Qnc that is released decreases as the o ffset distance r n ff set of the 
ignition region increases. This behavior is similar to that found by ITownsley et al.l (120071 ) in 
their 2-D simulations of the GCD model. 

The unburnt surface material in the initial collision region reaches T > 3 x 10 9 K 
but densities of only p ~ 1 5 — 10 6 g cm -3 , which are insufficient to produce a detonation 
(INiemeyer &: Woosleylll997l ). The collision, however, produces inward- and outward-directed 
jets. The outward jet ejects material at velocities t>j et ~ 40,000 km s _1 . The inward jet 
impacts the stellar surface, stalls, and spreads a little. This sequence of events compresses 
the hot (T > 3 x 10 9 K) material ahead of the jet to densities p > 1 x 10 7 g cm -3 , and in 
some cases even p > 2 x 10 7 g cm" 3 . These conditions exceed conservative conditions for 



initiation of a detonation (INiemeyer fc Woosleylll997l ; lRopke. Woosley. fc Hillebrandtll2007l ). 
As we will report in a later paper, subsequent 3-D simulations of the GCD model we have 
carried out that include nuclear burning outside the flame detonate, and inclusion of such 
burning does not significantly alter the time at which the detonation occurs compared to 
the first moment at which the simulations reported in the present paper reach conditions for 
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detonation ( Jordan et al.l 120081 ) 



The 3-D simulati ons of the GCD model show that, as in our previous 2-D simulations 
( ITownsley et al.l 120071 ). it is the kinetic energy originating from the breakout of the bubble 
of hot ash, imparted to the unburnt surface layers of the star by the inwardly-moving jet 
generated by collision of the surface flows that causes the unburnt material to achieve the 
conditions for detonation. This is illustrated in Figure [5j which shows the temperature and 
density in the collision region in the 18b42o6r simulation just before and just after detonation 
conditions are reached, and in Figure EH which is identical to the image in Figure 12(d), except 
that the star is rotated to show the very hot (T > 3 x 10 9 K) region at the head of the inward 
jet that has reached densities greater than 2 x 10 7 g cm -3 . 

We find that the unburnt material in the surface layers of the star reaches temperatures 
T > 2 x 10 9 K and densities p > 1 x 10 7 g cm -3 in all seven 3-D simulations we performed, as 
illustrated in Figure [6j The values of T max and p max closely match those in our 2-D cylindrical 
simulations for the same resolution and initial conditions (ITownsley et al.l 120071 ) . The small 
difference between the values of T max and p max in the 2-D and 3-D simulations for r oSset = 20 
km are within the uncertainties we expect in the calculations. Thus the results of our 2-D 
cylindrical simulations are a good guide to the results of our 3-D simulations for the range 
of initial bubble radii and offset distances, and the resolutions, that we have explored so far. 

In order to test the robustness of the GCD mechanism, we ran additional simulations in 
which we coarsened the resolution in the truncated cone encompassing the collision region 
from 8 km to 16 km, 32 km, and 64 km for offset distances of 40 km and 100 km. In all 
cases, the simulations reached the above conservative conditions for detonation. We conclude 
that, for the initial conditions investigated in this paper, the GCD model robustly achieves 
temperatures and densities necessary for detonation. 



4. Discussion 



4.1. Observational Properties of the 3-D GCD Model 



We have carried out 3-D simulations of the GCD model at a finest resolution of 6 km 
for initial offset distances of 42 km and 100 km, and a finest resolution of 8 km for initial 
offset distances of 20 km, 42 km, and 100 km. We find t hat t hese s imulations robustly reach 
the conditions necessary for detonation. IVladimiroval (120071 ) and iFisher et al.l (120081 ) have 
shown that the evolution of the bubble is self-similar - i.e., its evolution is independent of 
its initial radius nibble - provided that rabble < A c , where rabble is the initial bubble radius 
and A c is the minimum wavelength for the unstable Ray leigh- Taylor growth of flame surface 
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perturbations (see §2). We find that the bubbles in the simulations in which rabble = 25 
and 50 km at an offset distance r g set = 100 km do not exhibit self-similar behavior, even 
at early times; i.e., their size and shape at later times differ from each other and from those 
in the simulation for which rabble — 16 km (see Figure H]). This is expected, since in these 
two simulations the radius of the initial bubble does not satisfy the condition rbubbie < A c 
required for self-similar behavior. These two simulations produced larger values of E nnc , but 
still reached conservative conditions for detonation. Finally, we find that our 3-D simulations 
exhibit a correlation between E nuc and initial offse t distance (see Figure HI), confirming the 
correlation seen in our 2-D cylindrical simulations (jTownsley et al. 2007 ). 



Table 1 lists the amount of nuclear energy E nuc released up to a fiducial time of 2 s (which 
is approximately the time at which a detonation would occur) for the seven simulations 
reported in this paper. Also listed in Table 1 are estimates of the masses of heavy elements 
(i.e., iron-peak elements) and intermediate mass elements that are expected to be produced 
by a subsequent detonation, assuming that the mass of iron-peak elements Mh ea vy is that at 
densities above 1.5 x 10 7 g cm -3 and the mass of intermediate mass elements Mi nter is the 
mass at densities below this. Finally, Table 1 lists estimates of the expected total energy of 
the explosion, £7 to tai = E nuc + E hind , where we estimate E nuc from the masses of heavy and 
intermediate mass elements that are produced in each simulation, and E^ ind is the binding 
energy of the initial WD model. These results suggest that all of the simulations can be 
expected to produce large amounts of 56 Ni, and therefore very bright and energetic Type 
la supernova explosions, and small amounts of intermediate mass elements. Thus these 
simulations can explain only the brightest and most energetic of the Type la supernovae 
that are observed. 

Simulation 50bl00o8r, in which the initial bubble has a radius rbubbie — 50 km and 
therefore become immediately subject to a strong Rayleigh- Taylor instability, may crudely 
mock up what happens if ignition occurs simultaneously at a cluster of points located off- 
center in the core of the star. The simulation released more E nnc than did the other two 
simulations with the same offset distance but with smaller initial bubble radii, and thus 
suggests a plausible way in which the GCD mechanism could produce more pre-expansion, 
and therefore much less nickel, yet detonate - i.e., one way in which the GCD mechanism 
might account for less luminous Type la supernovae. However, this is an open question, 
which we plan to explore in a future paper. 

An essential aspect of the GCD model is that, while the nuclear energy released during 
the deflagration phase causes the star to expand prior to the detonatio n, it leaves the majority 



of the star unburn t and undisturbed, as w e have shown in earlier work (jPlewa. Calder. fc Lamb 



2004 ; iPlewa! 120071 ; iTownsley et al.ll2007l ) and in this work. The subsequent detonation phase 
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can therefore be expected to produce a smooth, stratified compositional structure in the 
interior of the star similar to that inferred from spectroscopic observations of Type la super- 
novae, and somethin g that 2 -D cylindrical and 3-D simulations of both t h e pure deflagration 
mode l jHoflich et all J2002h : Eeonard et all (b()05h : feadenes et all fcimh : IWang et al.l J2OO6 



20071 )]. and the deflagration to detonation (DDT) model [see, e.g. iGerardy et al.l (120071 )] 



have difficulty in doing. As is evident in Figures [T] and [2j the GCD model also produces tur- 
bulence and compositional inhomogeneities in the outermost layers of the star, which appear 
capable of matching pro perties inferred from observations of li ne polarization in th e optical 
f lWang et al Jbood . l2007h and line profiles in the NIR and MIR JGerardv et al.ll2007( fl Thus, 
the pure deflagration and DDT models predict an inhomogeneous, mixed composition in the 
core and a uniform composition in the outermost layers of the star - a compositional struc- 
ture that is opposite to that inferred from observations, while the GCD mechanism predicts 
a smoothly-stratified composition in the core and an inhomogeneous, mixed composition in 
the outermost layers of the star, which agrees qualitatively with the compositional structure 
inferred from observations. 



4.2. Comparison With Other Work 



Ropke. Woosley. fc Hillebrandtl (120071 ) have recently conducted an extensive set of 2- 
D cylindrical simulations and a few 3-D simulations of the deflagration phase of the GCD 
model. They find that the conditions for detonation are reached for a number of their 2-D 
cylindrical simulations for a variety of initial conditions. However, the 3-D simulations they 
performed did not reach conditions for detonation, whereas our 3-D simulations do. 

In an effort to understand the origin of this difference, we carried out 6-km and 8-km 
resolution simulations for exactly the same initial conditions as were used for one of the two 



Ropke. Woosley. fc Hillebrandtl (120071 ) 3-D simulations in which ignition was posited to occur 
at a single point: an initial spherical bubble of radius 25 km offset a distance of 100 km from 
the center of the star (see above) @ We also carried out an 8-km resolution simulation with 
an initial bubble radius rabble = 50 km and an offset distance r ff se t = 100 km. In all cases, 
the simulations reached conservative conditions for detonation, as we have described above. 



■"■It should be noted that one of the two events discussed in (jGerardv et al.ll2007l ) is a subluminous Type 
la supernova, whereas the GCD simulations reported in this paper produce very bright Type la supernovae. 



2 We did not simulate the other initial conditions for which iRopke. Woosley. fc Hillebrandtl (|2007T) did a 
3-D simulation positing a single ignition point (i.e., an initial bubble radius of 25 km and an offset distance of 
200 km) becau se these initial cond itions lie far above the rabble = K curve, and are therefore not physically 
self-consistent ( Fisher et al. 20081 ). 
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The results provide evidence of the ability of the GCD model to produce the conditions 
for detonation for a range of in i tial c onditions, but leave unanswered the question of why 
Ropke. Woosley. fc Hillebrandtl (120071 ) find that the criteria for detonation are reached for 
a range of initial conditions in their 2-D simulations but not in the 3-D simulations that 
they performed, whereas we find that th e criteria for detonati on are satisfied for a range of 
initial conditions in our 2-d simulations (ITownsley et al.l 120071 ) and in our 3-D simulations, 
as reported in this paper. 

We s how the values of r max and E nuc for both the current 3- D models and our previous 2- 



D mod els (ITownsley et aU20071 ). and the 2-D and 3-D models of lRopke. Woosley. fc Hillebrandt 
( 120071 ) in Figure [TJ This figure shows that there is a relation between T max in the collision 
region and E nuc . Such a relation is expected in the GCD model to the degree that larger 
values of E nnc produce more pre-expansion of the WD, and therefore less kinetic energy in 
the flow of hot bubble material over the stellar surface, leading to lower values of the tem- 
perature in the collision region. The amount of pre-expansion of the star can be expected 
to depend on when the nuclear energy is released, as well as how much is released. The fact 
that the relation between T max in the collision region and E nnc shown in Figured is relatively 
narrow suggests that the amount of nuclear energy that is released is the dominant factor in 
determining T max . 

The results of our 2-D cylindrical and 3-D simulations for initial conditions consisting 
of a single, small, spherical bubble offset a range of distances from the center of the star 
agree with each other, as previously noted, and lie on the relation between T max and E QVLC . 
So do the results of Ropke et al.'s 2-D cylindrical simulations for initial conditions consisting 
of a cluster of bubbles and for two tear-drop-shaped ignition regions located on opposite 
sides of the center of the star, as well as a single spherical bubble, all offset a range of 
distances from the center of the star. The results of their 3-D simulations for initial conditions 
consisting of a single bubble also follow the relation between T max and E nuc , but those for 
initial conditions consisting of a cluster of bubbles and for two tear-drop-shaped ignition 
regions located on opposite sides of the center of the star do not. In particular, their 3-D 
simulations starting with a cluster of bubbles release a low enough E QVLC that they should reach 
detonation conditions if the relation between T max and E nuc were followed. Most importantly, 



Ropke. Woosley. fc Hillebrandtl (120071 ) make no mention of the outward-directed jet, and 
especially the inward-directed jet, which we find plays a crucial role in achieving detonation 
conditions. Consequently, it is difficult to make direct comparisons between our simulations 
and theirs. 

The main conclusion we draw from Figure [7] is that the amount of E nnc released in 
Ropke et al.'s 2-D simulations is more, and in Ropke et al.'s 3-D simulations is much more, 
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than is released in our 2-D lTownsley et al.l (120071 ) and 3-D simulations, which release similar 



amounts of E nuc for the same initial conditions and resolution. Consequently, the kinetic 
energy of the surface flow in their 3-D simulations - as measured by T max - is much smaller 
and the simulations do not achieve detonation conditions. 

Without knowledge of the details of their simulations, it is difficult to know why the 
results of their 3-D simulations differ from those of their 2-D simulations, which - while 
releasing more E nvLC - lead to a surface flow similar to what we see in our 2-D and 3-D 
simulations. A likely reason is the different treatments of buoyancy-driven turbulent nuclear 
burning in their simulations and in ours. Our treatment assumes that the rate of buoyancy- 
driven turbulent nuclear burning depends mostly on the behavior of the flow at larger scales. 
As we have seen, this provides an explanation for why the results of our 2-D and 3-D 
simulations agree. It may also provide an explanation for why the results of Ropke et al.'s 
2-D simulations agree with our 2-D (and therefore also our 3-D) simulations: their treatment 
of buoyancy-driven turbulent nuclear burning uses the properties of the flow at scales above 
the grid scale of the simulation to determine the turbulent energy at subgrid scales, and 
therefore the increase in the nuclear burning rate (which is parameterized as an increase in 
the value of the flame speed) due to this turbulence. The fact that turbulence in 3-D leads 
to a cascade of smaller and smaller eddies, while turbulence in 2-D does not, means that 
their treatment of buoyancy-driven turbulent nuclear burning does not increase the rate of 
nuclear burning in 2-D, whereas in 3-D it will. The origin of the similarities between our 
2-D and 3-D results and Ropke et al.'s 2-D results, and the difference between these results 
and Ropke et al.'s 3-D results, are thus most likely due to differences in the treatment 
of buoyancy- driven turbulent nuclear burning. As we have noted above, the appropriate 
treatment of such burning is an open question, and is - as the above differences emphasize 
- an important topic for future study. 



5. Conclusions 

We have conducted a series of 3-D simulations of the GCD mechanism for several offset 
distances and resolutions. Conservative conditions necessary for detonation are robustly 
achieved in all cases. The initial conditions explored in this paper lead to conditions at 
detonation that can be expected to produce large amounts of 56 Ni and small amounts of 
intermediate mass elements. These particular simulations are therefore relevant only to high 
luminosity Type la supernovae. We find a correlation between the central density of the 
star at detonation and both the offset distance and the radius of the initial bubble. These 
correlations offer a possible explanation for the observed variation in nickel mass in Type 
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la supernovae. Finally, the uniform, homogeneous cores and the turbulent, heterogeneous 
composition of the outer layers of the stars at the time when the conditions for detonation are 
reached match the properties inferred from recent polarization, NIR, and MIR observations 
of Type la supernovae. 
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Table 1. Properties of 3-D GCD Simulations 



Label 


^bubble 


f offset 


Resolution 


-^miCjdef 


-^"heavy 


-Winter 


^total 


^det 




(km) 


(km) 


(km) 


(10 49 erg) 


(M e ) 


(M Q ) 


(10 51 erg) 


(s) 


16b20o8r 


16 


20 


8 


10.5 


1.00 


0.36 


1.50 


2.89 


18b42o6r 


18 


42 


6 


6.7 


1.18 


0.19 


1.57 


2.30 


16b40o8r 


16 


40 


8 


6.1 


1.20 


0.16 


1.58 


2.38 


16bl00o8r 


16 


100 


8 


3.2 


1.26 


0.10 


1.60 


2.02 


25bl00o6r 


25 


100 


6 


3.0 


1.27 


0.17 


1.60 


1.84 


25bl00o8r 


25 


100 


8 


3.1 


1.26 


0.11 


1.60 


2.01 


50bl00o8r 


50 


100 


8 


6.5 


1.11 


0.25 


1.55 


2.45 



Note. — This table gives the properties of the seven 3-D simulations of GCD models reported 
in this paper. 
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Fig. 1. — Images showing the hot ash and the star at different times for the 6- km resolution 
simulation of the GCD model starting from initial conditions in which a 25-km radius initial 
bubble offset 100 km from the center of the star. The images show volume renderings of 
the surface of the star (defined as the region in which p = (1.5 — 2.0) x 10 7 g cm -3 and the 
flame surface (defined as the surface where the flame progress variable <p\ = 0.5) at (a) 0.5 
s, soon after the bubble becomes Rayleigh- Taylor unstable and develops into a mushroom 
shape, (b) 1.0 s, as the bubble breaks through the surface of the star, and (c) 1.5 s, when 
the hot ash is flowing over the surface of the star, and (d) 1.7 s, shortly before the hot ash 
from the bubble collides at the opposite point on the surface of the star 
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Fig. 2. — Images showing very hot matter and the star at different times for the same 
simulation as in Figure [TJ The images are volume renderings of the surface of the star 
(defined as the region in which p = (1.5 — 2.0) x 10 7 g cm -3 and the regions where the 
temperature is very high (i.e., where T > 1.5 x 10 9 K) at (a) 0.8 s, when the bubble has 
become Rayleigh- Taylor unstable and developed into a mushroom shape; (b) 1.0 s, as the 
bubble breaks through the surface of the star; (c) 1.7 s, shortly before the hot ash from the 
bubble collides at the opposite point on the surface of the star, and (d) 1.84 s, the moment 
when the inward jet has compressed and heated stellar material ahead of it to detonation 
conditions. 
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Fig. 3. — Same image as that in Figure 2(d), except we have rotated the star to show the 
face of the inwardly directed jet where conditions for detonation are robustly achieved. 
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Fig. 4. — Nuclear energy released (E nnc ) as a function of time. Left panel: Comparison of 
6-km and 8-km resolution simulations for two offset distances. Note the close agreement 
between the two resolutions for both cases. Right panel: Comparison of 8-km resolution 
simulations for an initial bubble radius of 16-km and three different initial offset distances 
(20, 40, and 100 km) and for two other initial bubble radii (25 and 100 km) and an initial 
offset distance of 100 km. Note that the curves for initial bubble radii of 25 km, and 
especially 50 km - initial conditions that violate the requirement for self-similarity (see text) 
- are displaced from the other curves, even at early times. 
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Fig. 5. — Close-up view of 2-D slices of the region near the "south pole" of the star. The 
slices show the inward-directed jet produced by the collision of unburnt material ahead of 
the hot ash from the bubble in the 25bl00o6r simulation just prior to when the density of the 
material in the hot, inward-directed jet produced by the collision has reached its maximum 
value. The color shows the temperature, ranging from 1 xlO 9 - 5 xlO 9 K from blue through 
red. The density is indicated by contours. The yellow contour represents a density of 5 x 10 5 
gm/cm 3 , green 1 x 10 6 gm/cm 3 , purple 5 x 10 6 gm/cm 3 , red 1 x 10 7 g cm -3 , and black 2 x 10 7 
gm / cm 3 . 
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Fig. 6. — Maximum temperature T max and associated density in the fully refined truncated 
cone around the "south pole" of the star as a function of time for the five 8km resolution 3-D 
simulations we performed. The material flowing over the surface of the star enters the lower 
hemisphere at ~ 1.5 s and collides at ~ 2 s, at which point an inward-directed jet forms. 
Subsequently, the hot (T > 3 x 10 9 K) material in the jet impacts the surface of the star and 
becomes compressed, reaching densities p > 1 x 10 7 g cm -3 in all five of the simulations. 
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Fig. 7. — Locations in the (-E^uc, T max ) -plane of our 2-D and 3-D simulations of the GCD 
model and of Ropke et al.'s 2-D and 3-D simulations of the same model. The four filled dia- 
monds are (from left to right) the locations of our 25bl00o6r and 25bl00o8r simulations (for 
which the diamonds almost completely overlap) and our 16b 40o8r and 18b42o6r simulations. 
Note the correlation between E QVLC and T max reported by (jRopke. Woosley. fc Hillebrandt 
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